Big Bounce and inhomogeneities 
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The dynamics of an inhomogeneous universe is studied with the methods of Loop Quantum 
Cosmology as an example of the quantization of vacuum cosmological spacetimes containing grav- 
itational waves [Gowdy spacetimes). The analysis performed at the effective level shows that: (i) 
The initial Big Bang singularity is replaced (as in the case of homogeneous cosmological models) by 
a Big Bounce, joining deterministically two large universes, [ii) the universe size at the bounce is at 
least of the same order of magnitude as that of the background homogeneous universe, (iii) for each 
gravitational wave mode, the difference in amplitude at very early and very late times has a van- 
ishing statistical average when the bounce dynamics is strongly dominated by the inhomogeneities, 
whereas this average is positive when the dynamics is in a near-vacuum regime, so that statistically 
the inhomogeneities are amplified. 

PACS numbers: 04.60.Pp, 98.80.Qc 



The cosmological models based on classical General 
Relativity (GR) predict that the Big Bang singularity 
is the true beginning of the universe (boundary of the 
spacetime). This prediction, however, is believed not to 
be physical since, when one describes the early universe, 
GR is applied beyond its domain of validity. The quan- 
tum effects which dominate in this epoch are expected to 
resolve the singularity. 

This issue was studied recently in the context of Loop 
Quantum Cosmology (LQC) [l| (see also 0]). The anal- 
ysis of a simple model of a homogeneous and isotropic 
universe revealed surprising results Q: a large classical 
expanding universe was preceded by a(n also large and 
classical) contracting one, which bounced (deterministi- 
cally) and started to expand once the energy density of its 
matter content reached the Planck scale. The robustness 
of this result in more general situations was confirmed 
later [1, including the case of anisotropic universes. 

However, up to now these studies were mostly re- 
stricted to homogeneous spacetimes. Although a pre- 
liminary analysis of the effects of inhomogeneities was 
presented in [6], it was not known whether the modifi- 
cations to the universe dynamics predicted by LQG sur- 
vive in the presence of inhomogeneities. The possibility 
to answer this question arose when an LQG quantiza- 
tion scheme was formulated [7] for a class of cosmolog- 
ical spacetimes known as Gowdy universes [1]. These 
spacetimes, while still symmetric (they admit two spatial 
symmetries) include inhomogeneities that can be inter- 
preted as (linearly polarized) gravitational waves, thus 
having local degrees of freedom. To describe them, a hy- 
brid quantization scheme was applied: first the geometry 
was represented as the Fourier modes of a gravitational 
field (the linearly polarized wave) propagating in a ho- 
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mogeneous (Bianchi I) spacetime, next this Bianchi ge- 
ometry was quantized using loop techniques, while for 
the gravitational wave modes standard Fock quantiza- 
tion methods were employed after introducing a suitable 
time-dependent scaling. 

This construction paved the way to analyze the quan- 
tum dynamics of this inhomogeneous system; nonetheless 
the (field-like) complexity of the quantum configuration 
space makes the investigation of the genuine quantum 
evolution extremely difficult. This forces one to resort 
to classical effective dynamics [To| - a classical the- 
ory which incorporates the main effects of spacetime dis- 
creteness and which has been shown to accurately mimic 
the quantum evolution in the cases where it has been 
tested on so far. Here, we apply this technique to ana- 
lyze the quantum dynamical behavior of the Gowdy uni- 
verse, answering in particular the following questions: (z) 
Does the Big Bounce persist in the presence of inhomo- 
geneities? (m) If the answer is in the affirmative, does 
the Big Bounce occur in similar conditions as in homo- 
geneous models? {iii) And how does the structure of the 
gravitational wave modes change through the bounce? 

Let us start by specifying in more detail the physi- 
cal system which we consider here. The Gowdy uni- 
verses are vacuum spacetimes with compact sections of 
constant time which, in spite of possessing considerable 
symmetry, still contain local degrees of freedom. Namely, 
these spacetimes possess two spatial isometrics (two com- 
muting spacelike Killing vectors). The Gowdy universes 
can be classified by their spatial topology Q. The best 
studied case, on which we will concentrate, is that with 
the topology of a three-torus. This family of space- 
times provides a generalization of the homogeneous and 
anisotropic Kasner solution (for spatially flat topology) 
to include inhomogeneities which depend only on one spa- 
tial coordinate [ll|. The spacetime inhomogeneities can 
be interpreted as gravitational waves propagating in an 
homogeneous background spacetime. We will consider 
exclusively the simplest case of linearly polarized waves, 
in which, after a suitable gauge fixing, all the inhomo- 
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geneities can be described by a single metric field and 
the metric adopts a diagonal form globally 0, [13 • That 
field can be expanded in Fourier series exploiting the pe- 
riodicity in the only spatial coordinate on which it de- 
pends (this coordinate is cyclic for the studied topology) . 
Strictly speaking, the inhomogeneities are determined by 
the nonzero modes of this decomposition. The rest of 
gravitational degrees of freedom of the system describe a 
homogeneous universe on which the gravitational waves 
propagate. Specifically, this is a Bianchi I spacetime with 
three-torus topology. The classical solutions of this cos- 
mological model are known in exact form, and generically 
present a Big Bang singularity [Tst - 

In order to quantize the system we follow the pre- 
scription used in where the geometry (homoge- 
neous) degrees of freedom were quantized via loop tech- 
niques, whereas for the "matter" ones (in this case the 
gravitational waves with a convenient field parametriza- 
tion) standard (Schroedinger-like) methods were em- 
ployed The classical metric of the Bianchi I back- 
ground in the adopted gauge is diagonal and determined 
by a triple of scale factors ai{t), where i = 1,2,3. In 
the formalism used for quantization, the phase space for 
this background is coordinatized by the SU (2) (Ashtekar) 
connection = cMf/(27r) and the densitized triads 
Ei = piSi/{4:n'^) [ll|. All the information about the 
system is encoded in the canonically conjugate variables 
{c^,Pi), where \pi\ = \e^^''\ajak/{8TT) {e^^^ is the com- 
pletely antisymmetric unit tensor). Among the con- 
straints that GR imposes on the system, only two are not 
automatically satisfied in the introduced gauge: the gen- 
erator of 5*^ translations on the inhomogeneous spatial 
direction (which affects exclusively the inhomogeneities 
described by the gravitational waves) and the spatial av- 
erage of the Hamiltonian constraint. The latter can be 
written (up to a global constant) as C = Cbi + Cjt (where 
Cbi and Cjr stand for the homogeneous background part 
and the "matter" correction encoding inhomogeneities, 
respectively). Specifically, the background part equals 
Cbi = /s d^xl <lei{E)\-^l'^e''^''ElElF,j ^ where ^ is the 
curvature (field strength) of the connection A^. 

The quantization methods parallel those of LQG p^ . 
In a first step, the constraint is ignored. The basic ob- 
jects promoted to operators are the integrals of A\ along 
straight lines (holonomies) and those of Ef along square 
surfaces (fluxes) . The resulting kinematical Hilbert space 

is a product 7^"^'" = 0^^^ L'^{^b, d/is) where Mb 

is a Bohr compactification of the real line and is a 
Hilbert space for the inhomogeneous degrees of freedom. 

In the next step, the constraint C is promoted to (and 
solved as) a quantum operator. For this, Cbi has to 
be expressed in terms of holonomies and fluxes. In par- 
ticular, the term -Fijc entering Cbi is approximated by 
holonomies along small rectangular loops. Since the limit 
of the loop shrinking to zero does not exist in LQC, the 
rectangular loop is fixed (following [H]) by the require- 
ment that its physical area equals the lowest nonzero 
eigenvalue of an area operator defined in LQG. 



In the matter part of C, the degrees of freedom cor- 
responding to gravitational wave modes (conveniently 
scaled by a time-dependent function) are represented via 
fields and momenta operators combined into creation and 



annihilation operators {a\ 



m — ±1,±2,..., with 



standard commutation relations [12|. Therefore they 
form the standard Fock space 7i^" = T. 

The final quantum constraint C is a difference opera- 
tor in all three coefficients pi. In principle, one can find 
the (generalized) states annihilated by it, thus identifying 
the physical Hilbert space T^p'^''. This is indeed done in 
7]. However, the representation of any state of physical 
interest is complicated to the extreme by the presence of 
an infinite number of degrees of freedom. To be able to 
extract interesting physics out of the system we appeal 
here to the so-called classical effective dynamics. 

To derive this effective description (see Q for details), 
one replaces the basic operators (holonomies and fluxes) 
in the Hamiltonian constraint by their expectation val- 
ues, thus building back a classical Hamiltonian. Nonethe- 
less, some aspects of the quantum theory are preserved 
by leaving the lengths of holonomies as determined by 
the minimal area requirement discussed above. The final 
form of the constraint (after proper densitization) reads 
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Here, i — \ denotes the direction in which there exists 
spatial dependence, G is Newton's constant and M = 
\/ 2'\/37r7^pi (7 is the Immirzi parameter and Ipi — \J Gh 
the Planck length). We note that the variables 02 and 
03 are constants of motion, a fact which will considerably 
simplify the discussion of the effective dynamics. 

Given this effective constraint, it is straightforward to 
derive the complete set of equations of motion, namely 
the Hamilton-Jacobi equations. By integration, one then 
obtains the time evolution of the system from any initial 
point in the phase space. 

At this stage, it is worth pointing out that the above 
procedure does not take into account state-dependent 
parameters (like dispersions and higher-order quantities) 
which might significantly affect the dynamics. However, 
comparison of this effective approximation against the 
full quantum dynamics has been carried out in simpler 
models (i.e., homogeneous cosmologies with both mass- 
less and massive scalar fields), showing that the correc- 
tions arising from such parameters are negligible for the 
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FIG. 1: A dynamical trajectory for the case mmax = 5 (black) 
compared against the classical trajectories to which it con- 
verges in the distant future and past (red). Here, the value 
of Hi at the initial point (pi = p2 = 180/pi) is 0.657fi. The 
Immirzi parameter is set to its standard value, 7 = 0.2375 . . . 



states of physical interest (semiclassical at late times). 
This fact has been confirmed analytically in the simplest 
cosmological models Although the validity of the 

adopted effective description has not been tested in the 
model considered here, the commented results strongly 
support its ability to reproduce the behavior of the quan- 
tum system quite accurately. 

We probed the dynamics of the effective system us- 
ing Monte-Carlo methods. First, a large population (c.a. 
10^ points) of initial points was selected randomly in the 
phase space. In particular, apart from the restriction 
to satisfy the S'^-translation constraint, the values of a^n 
were generated randomly with Gaussian probability dis- 
tribution (centered at the origin). Owing to technical 
limitations, simulations were restricted to a finite number 
of nonvanishing modes (specifically, the cases considered 
were m < mmax = 5, 10, 20). To find the dynamical tra- 
jectories, the initial value problem consisting of the full 
set of equations of motion plus the chosen initial data was 
integrated via built-in adaptive methods of Mathematica. 
The analysis of the evolution was focused on two issues: 
(a) the existence of the bounce, and (b) the changes in 
the structure of the inhomogeneities through the bounce. 

In order to address the first of these issues, the evolu- 
tion of the inhomogeneous universe was compared against 
the evolution of its homogeneous counterpart, i.e. the 
universe determined by the same initial data for the ho- 
mogeneous degrees of freedom and all inhomogeneities 
set to vanish. The dynamical trajectories obtained nu- 
merically confirm the presence of the bounce in all three 
spatial directions (see Fig. [T]). Thus, the qualitative evo- 
lution picture stays the same as in homogeneous mod- 
els. Besides, analytical studies show that the values of 
P2 and P3 when they bounce do not depend on the a^'s; 
hence the bounces in the homogeneous directions occur 
at the same values of pi as for homogeneous universes. Fi- 
nally, to check the behavior of the bounce in the inhomo- 



FIG. 2: The set of points satisfying relation ([2)) for the case 
82 = 3750;pi, 63 = 2500;pi, and H^^JHI = 2 ■ 10"''. They 
form reflective boundaries for the dynamical trajectories. 



geneous direction, more extensive analytical/numerical 
studies were performed. Exploiting the properties of the 
equations of motion, the investigation was carried out 
separately (and using different methods) in two domains: 
6263 < and 6263 > 0. 

In the first case, it follows straightforwardly from ^ 
that the presence of inhomogeneities can only push the 
bounce away, that is, the bounce happens at a larger 
value of pi. Indeed, recalling that pi = {pi,C} and taking 
into account that pi = at the bounce, one gets the 
constraint satisfied at that point. 
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where _F is a positive definite function such that 
02,03; 0,0) = 0. The form of ^ implies imme- 
diately that, when O2O3 < 0, the contribution of F in- 
creases the value of pi (provided that ©2 and O3 remain 
the same). As a consequence, the bounce in the inho- 
mogeneous case always happens at larger universe sizes 
than in the homogeneous scenario. 

The case 0203 > required a detailed numerical anal- 
ysis, since now the two terms in brackets in the numer- 
ator of (HI) have opposite signs. The behavior of the 
points satisfying ([2]) with respect to the magnitude of 
the inhomogeneities is shown in Fig. [51 A feature which 
is worth noting is the existence of the throat which al- 
lows, in principle, that the dynamical trajectory may go 
"down it" through an infinite sequence of bounces and 
recollapses, reaching the singularity. However this would 
be a critical trajectory (the set of initial data leading to 
such evolution has zero measure in the phase space), thus 
generically the universe will bounce at finite pi. Besides, 
except perhaps for a very small subset of initial condi- 
tions (near the critical case), the bounce happens at a 
value of pi above certain bound, which is of similar order 
to the value found in the absence of inhomogeneities, p°. 
The bound provided by our analysis is « 0.13p°- 

Our study of the inhomogeneities was focused on dis- 
cussing how their energy distribution changes through 
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the bounce. This information is encoded in |am|(i)- 
Since, for large pi, these quantities converge to well de- 
fined limits, it is particularly interesting to investigate the 
corresponding amount of asymptotic change for large uni- 
verses, A\am\ = hmt^+oo |am| -limi^_oo \am\- Actually, 
once the full set of initial data is specified, the evolution 
of the universe (including the inhomogeneities) is deter- 
ministic, so that the value of A|am| is fixed. Nonetheless, 
A I am I acquires a stochastic nature if one restricts its at- 
tention to energies instead of amplitudes for each of the 
gravitational wave modes, thus ignoring the initial phases 
of the a„i's. Within the space of possible trajectories, on 
the other hand, one can distinguish two differentiated 
regimes: (i) inhomogeneity dominated, for which the dy- 
namics around the bounce is dominated by the content 
of gravitational waves, and (ii) near-vacuum, for which 
those waves introduce only small corrections to the vac- 
uum Bianchi I dynamics around the bounce. In case 
(i) our numerical analysis of a large population of uni- 
verses shows that, gcnerically, A|am| is an antisymmetric 
function of the initial phase of a^. Owing to this anti- 
symmetry, the expectation value of A|am| (namely, its 
average over the dependence on initial phases) is ensured 
to vanish. In case (ii), a similar analysis shows that the 
antisymmetric behavior of A|am| is generically lost and 
its average becomes strictly positive. Therefore, in the 
near-vacuum case the quantum geometry effects around 
the bounce pump energy into the inhomogeneities. 

It is worth noticing that, although our numerical anal- 
ysis has been performed in all cases for a finite number 
of nonvanishing modes rrimax (which plays the role of an 
UV cutoff on the gravitational waves) the results pre- 
sented here do not change when mmax increases. As a 
consequence, they remain valid when the full Fock space 
of inhomogeneities is considered. 

The results explained above constitute the first system- 
atic analysis of the dynamics of inhomogeneous space- 
times in LQC. In particular, they provide further sup- 
port to the bounce picture found earlier in homogeneous 



scenarios. The study reported here, and the methodol- 
ogy developed for it, paves the way for future analyses 
of the effective physics derived from LQG/LQC in more 
general spacetimes, opening an avenue, e.g., for the dis- 
cussion of the effects of quantum geometry in the process 
of gravitational collapse and black hole formation. 

Let us conchide clarifying that the presented analysis 
should not be treated as final. Firstly, our inhomo- 
geneous system has been investigated by means of an 
effective theory which does not take into account many 
quantum effects. To confirm the reliability of the results, 
one ought to repeat the analysis in the full quantum 
setting. Secondly, the method of construction of the 
minimal area loop used to define the quantum constraint 
C is not unique. There exist several prescriptions for the 
construction, giving different quantitative predictions of 
the dynamics. Therefore, for robustness, an investigation 
of the system constriictcd with those other prescriptions 
is also desirable. Finally, one has to remember that, in 
the studied system, there is only one inhomogeneous 
direction; therefore the results regarding change of 
inhomogeneities through the bounce cannot be directly 
applied to the kind of models considered in observational 
cosmology. To obtain reliable results verifiable against 
observations, one has to extend the analysis to models 
which have the energy level degeneracy characteristic of 
a spherical harmonics decomposition. This is the task 
that we plan to accomplish in future work. 

Acknowledgments: The authors are grateful to J. 
Cortez, L.J. Garay, J.M. Martin-Garcfa, and especially 
M. Martin-Benito. This work was supported by the 
Spanish MICINN Project FIS2008-06078-C03-03 and 
the Consolider-Ingenio 2010 Program CPAN (CSD2007- 
00042). D.B. acknowledges financial aid by the FPI Pro- 
gram of Madrid Regional Government, and T.P. by the 
ISP Program of CSIC (together with the ESF) and the 
Foundation for Polish Science grant Master. 



[1] M. Bojowald, Living Rev. Rol. 11, 4 (2008); A. Ashtekar, 

Nuovo. Cim. 122B, 135 (2007). 
[2] T. Cailleteau, A. Cardoso, K. Vandersloot, and D. 

Wands, Phys. Rev. Lett 101, 251302 (2008). 
[3] A. Ashtekar, T. Pawlowski, and P. Singh, Phys. Rev. 

Lett. 96, 141301 (2006); Phys. Rev. D 74, 084003 (2006). 
[4] A. Corichi and P. Singh, Phys. Rev. Lett. 100, 161302 

(2008). 

[5] D.W. Chiou, Phys. Rev. D 75, 024029 (2007); 76, 124037 
(2007). 

[6] M. Bojowald, M. Kagan, P. Singh, H. H. Hernandez, and 
A. Skirzcwski, Phys. Rev. Lett. 98, 031301 (2007). 

[7] M. Martm-Bcnito, L.J. Garay, and G.A. Mena Marugan, 
Phys. Rev. D 78, 083516 (2008). 

[8] R.H. Gowdy, Ann. Phys. 83, 203 (1974). 



[9] P. Singh and K. Vandersloot, Phys. Rev. D 72, 084004 
(2005). 

[10] V. Taveras, Phys. Rev. D 78, 064072 (2008). 

[11] B.K. Berger, Ann. Phys. 83, 458 (1974); 156, 155 (1984). 

[12] A. Corichi, J. Cortez, and G.A. Mena Marugan, Phys. 

Rev. D 73, 084020 (2006). 
[13] V. Moncricf, Pliys. Rev. D 23, 312 (1981). 
[14] A. Coriclii, J. Cortez, and G.A. Mena Marugan, Phys. 

Rev. D 73, 041502(R) (2006); A. Corichi, J. Cortez, G.A. 

Mena Marugan, and J.M. Velhinho, Classical Quantum 

Gravity 23, 6301 (2006). 
[15] M. Martm-Bcnito, G.A. Mena Marugan, and T. 

Pawlowski, Phys. Rev. D 78, 064008 (2008). 
[16] C. Rovelli, Living Rev. Rel. 11, 5 (2008). 



